Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 
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Abstract. We report the detection of a complex of narrow X-ray absorption features from the low-mass X-ray 
binary GX13+1 during 3 XMM-Newton observations in 2000 March and April. The features are consistent with 
being due to resonant scattering of the Ka and K/3 lines of He- and H-like iron (Fe xxv and Fe xxvi) and H-like 
calcium (Caxx) Ka. Only the Fexxvi Ka line has been previously observed from GX 13+1. Due to the closeness 
in energy the Fexxv and Fexxvi Kf3 features may also be ascribed to Nixxvn and Nixxvin Ka, respectively. 
We also find evidence for the presence of a deep (r ~ 0.2) Fe xxv absorption edge at 8.83 keV. The fits also 
require the presence of a broad emission feature whose energy and width are poorly determined, partly due to 
the presence of the deep Fe Ka features which severely cut into the feature and partly due to fit differences when 
using different XMM-Newton instruments. The equivalent widths of the lines do not show any obvious variation 
on a timescale of a few days suggesting that the absorbing material is a stable feature of the system and present 
during a range of orbital phases. 
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1. Introduction 



persistent low-mass X-ray binary 
1 (4U 1811-17) is an X-ray burst 
Matsuba 



1995) with radio 



The bright and 
(LMXRB) GX13 
source (Fleischman 1985 
(Grindlay & Seaquist 1986) and infrared counterparts 
(Naylor et al. |l99l| ; Charles & Naylor |1992| ). It displays 
a 15% modulation in X-ray intensity with a period of 
24.7±1 days. The X-ray intensity is anti-correlated with 



the hardness ratio (Corbet 1996). Recent estimates of the 
spectral type of the companion (an evolved late- type K5 in 
star; Bandyopadhyay et al. 1999| ), together with the bi- 
nary system parameters are consistent with this modula- 
tion being the orbital period, although further confirma- 



tion is required (Bandyopadhyay et al. 2002). GX13+1 
has been classified as an atoll source, but in some respects 
its properties are m ore similar to those of the Z-sources 
(Hom an et al. 1998| ; Schulz et al. 198S ; Bandyopadhyay 



et al. 1999 ), having persistent radio emi ssion at a similar 
level to the Z-sources (Fender & Hendry ^000| ) , and a high 
accretion rate of ~0.1 Eddington. A quasi-periodic oscil- 
lation (QPO) was discovered at 57-69 Hz by Homan et al. 
( |l998| ). Distance estimates for GX 13+1 range from 4.7 to 
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7.0 kpc (Christian & Swank 1997). Infrared spectroscopy 



favors a distance of 7+1 kpc (Bandyopadhyay et al. 1999) 



A spectral survey of Fe-K emission from 20 LMXRB 
with ASCA is reported in Asai et al. ( 2000 ). Significant 
iron emission was detected from GX 13+1 with a line cen- 
troid of 6.44 ± 0.05 keV, an equivalent width (EW) of 
37 ±1° eV and a line width cr~130±40 eV. The X-ray con- 
tinuum could be modeled with a two-component model, 
consisting of a multicolor disk-blackbody with a tempera- 
ture, kT, at the inner disk radius of 0.76 keV and a black- 
body with a kT of 1.34 keV, absorbed by a column density, 
N K , of 2.9 10 22 atom cm -2 . The source flux (1-10 keV, cor- 
rected for absorption) was 1.2 10 -8 erg cm~ 2 s^ 1 , which 
translates into a luminosity of 710 37 erg s _1 (at 7 kpc). 



Recently, Ueda et al. ( gOOlJ ) re-analyzed the ASCA Solid- 
State Imaging Spectrometer (SIS) observation of GX 13+1 
and discovered that the spectral region around Fe is more 
complex than had been realized. Their fits required the 
presence of a narrow (a < 70 eV) absorption line at 
7.01 ± 0.03 keV with an EW of 35+8 eV, which they 
interpreted as resonant scattering of the Ka line from 
Fexxvi ions with a column density <;10 18 cm -2 , a narrow 
(a < 220 eV) emission line at 6.42 + 0.08 keV with an EW 
of 19 ± 8 eV and an absorption edge at 7.61 ± 0.13 keV 
with an optical depth of 0.13 ± 0.05. This was the first 
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Fig. 1. The EPIC PN spectra for the 3 on-axis observations of GX13+1. In the left panels, the best-fit continuum 
model consisting of absorbed multi-color blackbody and blackbody components, together with a broad emission feature 
at 6.8 keV is shown with the normalizations of the narrow discrete absorption features, and the optical depth of the 
absorption edge, set to zero. On the right, expanded views of the residuals (in units of standard deviation) in the 
5.5-9.5 keV energy range. The four absorption lines, a broad emission line at 6.8 keV which is strongly affected by the 
two strongest absorption features and the absorption edge at 8.83 keV are clearly visible 



detection of an Fe absorption line in an X-ray binary 
known to contain a neutron star. Previously, the only X- 
ray binaries known to exhibit such Fe absorption features 
were two superluminal jet sources GRO J1655-40 (Ueda et 
al. |1998|; Ya maoka et al |200l|) and GRS 1915+ 105 (K otani 



et al. EOOOt Lee et al. pOOll ). As Ueda et al. (|200l|) com 



ment, the interpretation of the edge is not simple and it 



could be a combination of an absorption edge from mod- 
erately ionized Fe together with line absorption from a 
number of prominent Ka or K/3 lines of Fe and Ni. 

XMM-Newton is significantly improving our under- 
standing of the LMXRB emission processes, accretion 
geometries and the physical conditions of the material 
through which the X-ray emission propagates. Reflection 
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Table 1. XMM-Newton on-axis observation log of GX 13+1. All observations used the medium filter 
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Rev. 


Obs. ID 


Inst. 


Start time (2000) 


Net Exp. 


Mode 






EPIC 


(dy mon hr:mn) 


(ks) 




0056 


0101 


MOS1 


30 Mar 15:13 


5.9 


Full Frame 






PN 


30 Mar 14:47 


2.5 


Timing 


0057 


0901 


MOS1 


01 Apr 06:28 


4.9 


Full Frame 






PN 


01 Apr 05:52 


2.1 


Timing 


0057 


1001 


MOS1 


01 Apr 09:54 


4.6 


Full Frame 






PN 


01 Apr 09:28 


2.3 


Timing 



Grating Spectrometer (RGS) observations of the ecl ipsing 
and dipping LMXRB EXO 0748-676 (Cottam et al. |200l|) 
revealed several discrete absorption and emission features 
from ionized neon (Nex Lya, Neix Lya), oxygen (Ovill 
Lya, Ovu He-like complex) and nitrogen (Nvn Lya). 
These line features are broadened (velocity widths of 
1000-3000 km s _1 ), but with no evidence for a velocity 
shift (<300 km s _1 ). No obvious dependence on orbital 
phase was found. 

During an XMM-Newton observation of the eclipsing 
and dipping LMXRB MXB 1659-298 resonant absorption 
features identified with O vin ls-2p, ls-3p and ls-4p, Ne X 
ls-2p, Fe XXV ls-2p, and Fe xxvi ls-2p transitions, to- 
gether with a broad Fe emission feature at 6.47lg keV 



were discovered (Sidoli et al. 2001). The line widths are 
unresolved in the RGS and correspond to velocities of 
<600 km s" 1 for O vin ls-2p and <2300 km s" 1 for 
the Fe xxvi Ka line observed using the European Photon 
Imaging Camera (EPIC). The absorption features occur 
through a wide range of orbital phases, so the absorb- 
ing material is most likely located above, or below, the 
accretion disk, in a cylindrical geometry with an axis per- 
pendicular to that of the disk. 

Similar ab sorpt ion features have been reported by 
Parmar et al. ( 2002 ) during an XMM-Newton observation 
of the LMXRB X 1624-490, sometimes called the "Big 
Dipper". Here the lines are identified with the resonant 
Ka absorption of Fexxv and Fexxvi, with Nixxvn Ka 
and Fe xxvi K/3 possibly detected. Again, no obvious de- 
pendence on orbital phase was observed, except during 
an absorption dip. Here, we report the discovery of dis- 
crete X-ray absorption features from highly ionized Fe 
and Ca (and possibly Ni) in XMM-Newton EPIC spectra 
of GX13+1. 

2. Observations 



The XMM-Newton Observatory (Jansen et al. |2001| ) in- 
cludes three 1500 cm 2 X-ray telescopes each with an EPIC 
at the focus. Two of the EPIC imaging spectrometers 
use MOS CCDs (Tur ner et al. |200l| ) and one uses a PN 
CCD (Striider et al. 2001 ). The region of sky contain- 
ing GX13+1 was observed by XMM-Newton a number 
of times as part of the calibration programme. Three on- 
axis observations were performed (see Table 1). Since an 
analysis of two off-axis observations (2000 April 1 13:34 



and April 3 13:34 UTC) give similar spectral results, we 
discuss only the better calibrated on-axis observations. In 
order to minimize the effects of pile-up the PN and MOS2 
were operated in their Timing Modes. The MOS1 CCD 
was operated in its normal Full Frame Mode and the data 
are strongly affected by pile-up (when two photons are 
detected closely in time and the energy is incorrectly as- 
signe d to be the sum of the individual energies; sec Striider 
et al. 2001). The effects of pile-up in MOS1 were mini- 



mized by extracting events in an annulus outside the 30" 
radius core of the GX13+1 point spread function (PSF). 
Since there are far fewer MOS events than PN ones, the 
MOS1 data are only used to verify the overall shape and 
energy calibration of the PN data. 

Raw data products were extracted from the public 
XMM-Newton archive and then reprocessed using version 
5.2 of the Science Analysis Software (SAS), before be- 
ing further filtered using xmmselect. Only X-ray events 
corresponding to patterns (single pixel events) were se- 
lected. Background spectra were obtained from source free 
regions taken from the same observations. For the Timing 
Mode observations, appropriate public response matrices 
are not yet available. Thus, the standard PN matrix for 
Full Frame Mode has been adopted. This slightly affects 
the modeling of the continuum, and possibly the energies 
of the lines (F. Haberl, priv. communication), but is un- 
likely to affect estimates of their EWs. We concentrate 
here on the analysis of the narrow absorption features ev- 
ident in the spectra and defer a detailed discussion of the 
continuum properties to a later paper. 

3. Results 

A total of 6 EPIC spectra of GX13+1 (3 MOS1 and 3 
PN) from the on-axis observations reported in Table 1 
were extracted. These were then rebinned to oversample 
the FWHM of the energy resolution by a factor 3 and to 
have additionally a minimum of 20 counts per bin to allow 
use of the x 2 statistic. In order to account for systematic 
effects a 2% uncertainty was added quadratically to each 
spectral bin. The photo-ele ctric absorption cross sections 
of Morrison & McCammon( |l983|) are used throughout. All 
spectral uncertainties are given at 90% confidence, unless 
indicated otherwise. 

The overall 2.0-10 keV continua were investigated us- 
ing the same model as applied to the ASCA spectrum re- 
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Table 2. Parameters of the absorption lines detected in the XMM-Newton spectra of GX 13+1. A continuum model 
consisting of absorbed blackbody and disk-blackbody components was used 



Component 


Parameter 


Obs. 


0101 


Obs. 0901 


Obs. 1001 






MOS1 


PN 


MOS1 


PN 


MOS1 


PN 


Caxx abs 


Sitae (keV) 


4 1 Cj +0-05 
4 - i0 -0.05 


4 15 +0-3 

4 ' 10 -0.05 


4 16 +003 

*±.J-U _o 03 


4-18 + _° f 5 


4 iq+0-04 

^•J-O _o 03 


4 1Q+004 
_o 05 


line 


cr (keV) 


<0.05 


<0.1 


<0.1 


<0.2 


<0.05 


<0.1 




EW (eV) 


_ 7 fi +7.i 

' -7.6 


q q +1-9 
O.O _ 2 i 


-H-3 + 7 


-4.3l 2 23 


-7.7 tH 


-6.0+ 2 2 6 9 


Fexxv Ka abs 


Sitae (keV) 


6 7S+ 04 

u - ' ° -0.04 


°- ' d -0.02 


6.76 +°o1 

— U.U5 


6-75 

— O.Oz 


6.72+°°* 

— 0.06 


6.72+°'°? 


line 


cr (keV) 


<0.1 


<0.1 


<0.1 


<0.1 


<0.05 


<0.04 




EW (eV) 


oq o +13.2 
^•J.z _ 19 2 


9q 9 +4.6 
zo - z -15.1 


-26.7+173 


37.5 yj o 


-17-2 


q 9 9 +10.1 

8 2 


Fe xxvi Ka abs 


E linc (keV) 


7 00 + 002 
' -UU -0.02 


7 00 +°' 03 
'• uu -0.10 


7.06+2-gt' 

— U.U4 


7.03+S-jg' 


7.03 +°°3 

— U.Uo 


7.05+°'°i 

— U.Uz 


line 


cr (keV) 


<0.05 


<0.04 


<0.05 


<0.05 


<0.05 


<0.06 




EW (eV) 


n o+lS.O 
jr.o _ 18 5 


-36.1 


-57.0 1 2 °| 


-56.5t« 4 5 6 


-39.6 tlH 


-39.4+^ 


Fexxv K/3 abs 


Eunc (keV) 


7.85 fixed 


7 gn+ 01 

1 .»U _ 04 


7.85 fixed 


7 Q7 +0.05 
-0.03 


7 S/ i +0.04 

7-84 _ 04 


7 or +0.05 
' - 00 -0.02 


feature 


cr (keV) 


<0.05 


<0.05 


<0.05 


<0.04 


<0.05 


<0.05 




JiW (eV J 


oo i +28.1 
-28.1 _ 30 g 


Ol A -f-5.1 
- 2L4 -6.7 


—64 + 33 - 8 
°*- U -35.6 


-29.41^ 


7 q q + 29.6 

' OiJ -30.7 


_9Q Q+3.6 


Fexxvi K/3 abs 


Eunc (keV) 


8.26 fixed 


8-20 ±g;S 


8.26 fixed 


8 25 + ' 05 

-0.03 


o- 11 -0.01 


o 9 q +0.07 

°- zo -0.03 


line 


cr (keV) 


<0.05 


<0.05 


<0.05 


<0.05 


<0.05 


<0.06 




EW (eV) 


>-41.1 


-25.0 j*° 


ir. 7 +15.7 
10 -' -56.6 


9 7 q +8.4 
^ ' -9.6 


-69.1 ± 3 £;l 


_1Q q+6-9 
1».0 _ 10 . 4 


Fexxv edge 


Sitae (keV) 


8.83 fixed 


8.83 fixed 


8.83 fixed 


8.83 fixed 


8.83 fixed 


8.83 fixed 




T 


<0.40 


19 +°° 2 

U - ±J -0.02 


<0.42 


U-J-^ -0.08 


<0.32 


n 90 +0.04 
U - 28 -0.02 


Fexxvi edge 


Sitae (keV) 


9.28 fixed 


9.28 fixed 


9.28 fixed 


9.28 fixed 


9.28 fixed 


9.28 fixed 




r 


<0.22 


<0.04 


<0.50 


<0.03 


<0.20 


<0.01 


Fe broad emission 


E line (keV) 


6 24 + ' 27 

D - z ^ -0.24 


6.80l°; 3 9 


6.80 fixed 


6 85 + 11 

0-80 _ 05 


6.80 fixed 


f> 77 +0.02 
°- ' ' -0.09 


line 


cr (keV) 


0.8 ±0.2 


0.37 1° 2 ? 


0.3 fixed 


S5+ 014 

U.OO_ 08 


0.3 fixed 


26+ 13 

U.ZD_ 04 




EW (eV) 


260 +27 ° 

ZDU -160 


91 + 9 

al -30 


<130 


961S 


1701?° 


141+30 



ported in Ueda et al. ( 2001 ) consisting of absorbed black- 
body and multicolor disk-blackbody (Mitsuda et al. 1 



components. Fits were made to all 6 spectra individually. 
Examination of the fit residuals shows that, especially in 
the 3 PN spectra, broad positive residuals are present near 
6.7 keV, requiring the addition of a broad Gaussian line in 
all 3 fits with EWs of ^100 eV and widths, cr, in the range 
^200-400 eV. The continuum parameters derived from 
all spectra are similar. For the multicolor disk-blackbody 
component kT is in the range 0.5-0.7 keV while the kT 
of the blackbody is 1.1—1.2 keV. Both components are ab- 
sorbed by an N H , of 3.1-3.8 1 22 atom cm" 2 . Without the 
addition of the emission and absorption features discussed 
below the fit quality is rather poor with a xt 01 3.5-4.6. 
The addition of emission and absorption lines (modeled 
by Gaussians) successively decreases the xl to 2.0-2.4 (see 
below). The spectral parameters for the multicolor disk- 
blackbody and blackbody are consistent with th ose ob - 
tained from the ASCA observation by Asai et al. ( [2000 ). 



Two strong absorption features are observed in all 6 
spectra at around 6.7 and 7.0 keV (see Table 2). These 
features are consistent with being due to resonant K shell 
absorption from highly ionized iron Fexxv and Fexxvi, 
respectively. The deeper 7.0 keV feature is almost cer- 
tainly the a bsorp tion feature observed from GX 13+1 by 
Ueda et al. (|200l| ) using ASCA. The addition of both fea- 
tures is highly significant. Including the 7.0 keV line re- 
duces the xl by 0.6-0.8 while the addition of the 6.7 keV 



feature reduces the xl by 0.3-0.6. Two other absorption 
features at around 7.8 and 8.2 keV are also detected here 
for the first time in GX 13+1 in all the PN spectra, while 
in the MOS1 spectra they are visible only in one obser- 
vation (Obs. ID 1001), consistent with the lower effective 
area of this instrument, compared to the PN. The ad- 
dition of these features reduces the xt by 0.15-0.3 and 
0.04-0.07 for the 7.8 and 8.2 keV features, respectively. In 
both the MOS1 and PN spectra an absorption feature at 
~4.1 keV is also evident. This is consistent with being due 
to resonant Ka absorption from H-like calcium (Caxx). 
All the detected absorption lines and their parameters are 
listed in Table |^. We note that there is a strong interplay 
between the parameters of the broad emission line and 
the EWs of the two strong absorption features since these 
features cut strongly into the emission line. This is most 
clearly seen in Fig. [l| which shows the best-fit PN spectra 
together with expanded plots of the spectral regions that 
include the Fe features. 

Line energies measured using the PN Timing Mode 
spectra were expected to be slightly higher than their 
"true" values due to uncertainties in the calibration of 
Timing Mode. However, a comparison of the energies de- 
rived from the MOS1 and PN spectra for the two strongest 
absorption lines (Fe xxvi and Fe xxv Ka) shows that they 
are consistent at 90% confidence. We therefore present 
in Table § the PN energies without any adjustment. The 
identification of the fainter absorption features at 7.8 and 
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8.2 keV, is uncertain. They may be tentatively identified 
with Fexxv (7.88 keV) and Fexxvi K/3 (8.26 keV), or 
with Ni xxvn (7.80 keV) and Ni XX VIII Ka (8.10 keV), 
respectively (e.g., Kotani et al. 200C). 

Close inspection of the residual plots (Fig. [I]) reveals 
the presence of an emission feature near 8.5 keV. A fit 
to this feature with a Gaussian model gives an energy of 
8.51 ± 0.05 keV and an EW in the range 25-35 eV. Since 
this energy is close to those of known PN background fea- 
tures due to Zn Ka and Cu K/3 (an even stronger line is 
present in the background spectrum near 8.05 keV from 
Cu Ka) we investigated whether this feature could origi- 
nate from the background. In the PN Timing Mode, data 
from the on-axis CCD are continually clocked across the 
target region and data are extracted with only one dimen- 
sion of spatial information. We searched all the publically 
available XMM-Newton observations but were unable to 
find any blank field Timing Mode data from which to 
extract a background spectrum. Since the internal back- 
ground is spatially very structured, particularly in the K- 
fluorescent lines of metal components of the camera body, 
we constructed a background spectrum from the identical 
region covered by Timing Mode from Imaging Mode data, 
taking care to apply the appropriate exposure time correc- 
tion. From this spectrum we find that the observed 8.5 keV 
feature is a factor 100 more intense than the background 
features. We therefore conclude that the observed feature 
is not a background artifact. The energy of the emission 
feature does not correspond to any prominent emission 
line from an abundant element. However, it is just below 
that of the He-like Fe edge at 8.83 keV. We therefore at- 
tempted to model this feature as an absorption edge. The 
inclusion of an edge at 8.83 keV accounts for the excess 
around 8.5 keV and decreases the xl to 1.6-1.7. The best- 
fit parameters for all the discrete features can be found 
in Table where the upper limits to the presence of a 
H-likc Fc edge (at 9.28 keV) are also reported. A narrow 
emi ssion feature near 6.4 keV was reported by Ueda et 
al. ( 2001 ). Addition of this narrow feature with an energy 
close to 6.4 keV resulted in a decrease of the \ 2 °f 4-7, 
which is not significant. 

We note that the parameters derived for the broad 
emission feature from the PN and MOS1 are inconsis- 
tent. The energy derived from the MOS is lower, while 
the width of the line is significantly broader. This is not 
caused by the fitting method (e.g. a local minimum) since 
the addition of the line with the PN-derived parameters 
to the MOS spectral model results in a clearly worse \ 2 
(Ax 2 > 100), and the difference between the model and the 
MOS data is up to 10-15%, i.e. larger than the calibration 
uncertainty. Since a large fraction of the line emission is 
removed by the narrow absorption components it is dif- 
ficult to reliably constrain the emission component. We 
therefore await an improved EPIC calibration before the 
parameters of the broad emission line can be accurately 
determined. 



4. Discussion 

There is a growing evidence that narrow absorption lines 
from highly ionized material are a common feature of 
LMXRB spectra. They were first detected from the su- 



perluminal jet sources GROJ1655— 40 (Ueda et al. 1998; 
Yamaoka et al. |200l|) and GRS 1915+105 (Kotani et 



al. |2000| ; Lee et al. p00l[ ) with ASCA. Absorption lines due 
to Fe XXV and Fe xxvi were observed from GRO J1655— 40 
and did not show any obvious dependence of their EWs 
on orbital phase. Both X-ray and optical observations in- 
dicate that GRO J 1655— 40 is most probably a binary sys- 
tem viewed at a high inclination angle of ~70° (Kuulkers 
et al. 1995; van der Hooft et al. 1997; Greene et al. 



pool] ) 



ASCA observations of GRS 1915+105 revealed, in 
addition, absorption features due to Caxx, Nixxvn and 
Ni xxviii. This source has been observed also with the 
Chandra HETGS which detected absorption edges of Fc, 
Si, Mg, and S as well as resonant absorption features 
fro m Fe xxv and Fexxvi and possibly Caxx (Lee et 
al. 2001). The recent discovery of similar narrow iron 
absorption features in X-ray binaries containing neutron 
stars (GX 13+1, Ueda et al. |200l| ; MXB 1 658- 298, Sidoli 



et al. 2001 and X 1624-490, Parmar et al. 2002), rules out 



the possibility that such features are peculiar to superlu- 
minal jet sources and related in some way to the jet forma- 
tion mechanism. Chandra HETGS observations have also 
convincingly detected narrow absorption features due to 
Ne ix and Fe xvin from the transient black hole candidate 



XTE J1650— 500 (Miller et al. [2002]). These lines are may 
be due to absorption in an ionized accretion disk atmo- 
sphere or wind. 

The XMM-Newton observations reported here reveal 
a much more complex range of emission and absorption 
components from GX 13+1 than even the ASCA analy- 
sis of Ueda et al. ( (2001 ) had found. In addition to the 
Fe xxvi Ka absorption feature additional narrow absorp- 
tion features which may be identified with Ca xx and 
Fe xxv Ka and with Fe xxv and Fe xxvi K/3 are evident 
in the XMM-Newton spectra. In addition a deep edge at 
8.83 keV is clearly required. We note that the feature at 
7.61 ± 0.13 keV which Ueda et al. fl200l| ) model as an 
edge is more clearly seen to be a line absorption feature 
in the XMM-Newton PN spectra (which we identify with 
Fe xxv K/3). It is interesting to speculate as to why a 
deep (t ~ 0.2) edge from Fe xxv is evident whereas the 
PN upper-limits to the optical depth for an Fe xxvi edge 
at 9.28 keV of ^0.04 exclude a similar feature. This is in 
contrast to the corresponding absorption lines, where the 
EW of the Fe xxvi transitions exceeds that of Fe xxv in 
all the spectra. 

In order to convert the observed EWs of the lines into 
ion column densities, we assumed t he cu rve of growth 
analysis performed by Kotani et al. ( 200C| ). In this way, 
the column densities of the lines can be estimated for a 
range of kinematic temperatures of the absorbing mate- 
rial, where this temperature includes contributions from 
thermal motions as well as any bulk motions or turbu- 
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Log kT (keV) 



Fig. 2. Fe column densities for the two Ka absorption 
lines from the first (ID. 0101) MOS1 (squares are for 
Fexxvi, triangles for Fexxv) and PN observations (dia- 
monds are for Fe xxvi, stars for Fe xxv) . For the other two 
observations (ID. 0901 and ID. 1001) similar results have 
been obtained. The horizontal line marks the upper limit 
to the Fe column density corresponding to a Thomson op- 
tical depth equal to unity. A small horizontal shift has 
been applied for clarity. Column densities are in units of 
cm -2 
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Fig. 3. Results of the curve of growth analysis for the 
Caxx line, using only results from the first observation 
(Obs. ID 0101). Similar values are obtained from the other 
two observations. Squares corresponds to PN results, dia- 
monds to MOS1. A small horizontal shift has been applied 
for clarity. Column densities are in units of cm~ 2 



lence. An Fe abundance of 4.7 10 5 was assumed (Anders 



& Grevesse 1989). The column densities for Fexxv and 
Fexxvi, calculated for a range of assumed temperatures 
(0.1, 1, 10, 100 and 1000 keV), are shown in Fig. | for 
only one observation, for clarity. Indeed, similar EWs have 
been measured in all the three observations. A lower kine- 
matic temperature requires a higher ion column density, 
and hence a higher Nh- A constraint on the absorbing 
Fe column density can be derived from the fact that the 



plasma should be optically thin to Thomson scattering 
(N H < 1.5 10 24 atom cm" 2 ), other wise the absorption 
lines would be strongly diminished. Assuming a solar Fe 
abundance, implies an Fe column density <710 19 cm~ 2 , 
which translates into a temperature ^10 keV. The up- 
per limits to the velocity widths of the two narrow lines of 
<1800 km s" 1 (Fexxv) and <1900 km s -1 (Fexxvi), pro- 
vide not very constraining upper limits to the kT of 600- 
650 keV. The EWs of the highly ionized Fe K lines (and 
thus the ions column densities) are similar to those mea- 
sured in the eclips ing, dipping LMXRB MXB 1658-298 
(Sidoli et al. |200l| ) with XMM-Newton, and a factor -3 
higher than that observed in the LMXRB X 1624-490 



(Parmar et al. 2002 ). The results of the curve of growth 
analysis for the Caxx line are shown in Fig. |^. The ion- 
ized calcium column density ranges from 10 16 cm~ 2 to 
10 20 cm" 2 , depending on the kinematic temperature. No 
useful constraints come from the upper limit to the veloc- 
ity width of <3700 km s _1 , that translates into a temper- 
ature of ^1900 keV. 

The line centroids of the Fe xxv and Fe xxvi Ka do not 
show evidence for any velocity shifts, except in the second 
observation (Obs. ID 0901) where blue-shifts of -2700 ± 
2200 km s" 1 (Fexxv Ka) and -3900 ± 1700 km s" 1 



(Fexxvi Ka) are evident. In this 
blue-shift in the Caxx feature of - 



same observation, a 
5100 ± 2200 km s" 1 
is also detected. For most observations, the MOS energy 
calibration is correct to within —5 eV, or a velocity shift 
of -220 km s" 1 at 7 keV and -400 km s" 1 at 4 keV. 
Individual observations have been analysed where the Mn 
Ka internal calibration line appears to be discrepant by 
—30 eV, but these are all associated with anomalous in- 
strument temperatures resulting from spacecraft irregu- 
larities. Such conditions did not pertain for the observa- 
tions reported here. The absorption lines at higher en- 
ergy cannot be unambiguously identified. These features 
are only clearly detected in the MOS1 in the third ob- 
servation (Obs. ID. 1001). The measured line energy of 
7.84±0.04 keV can be associated with both the Fexxv K/3 
(at 7.88 keV) and the Nixxvn Ka (at 7.80 keV) feature. 
The absorption line at —8.1 keV can be similarly associ- 
ated with either an Fexxvi K/3 transition (at 8.212 keV) 
or with a Nixxvill Ka feature at 8.034 keV. 

Since the orbital period of GX 13+1 has not been con- 
firmed, we cannot reliably estimate a geometry and size 
for the absorbing plasma around the central X-ray source. 
If the 25 day period is confirmed to be orbital, we note 
that the XMM-Newton observations reported here only 
cover a small fraction (—10%) of an orbital cycle. All we 
can say is that all the narrow absorption lines are consis- 
tent with being present with the same EW in all analyzed 
spectra. This implies that the absorbing plasma is stable 
in the line of sight over timescales of few days suggest- 
ing a stable geometry, as in the case of MXB 1658—298 
and X 1624-490, and excluding the possibility of absorbing 
clouds temporarily passing in front of the central source 
producing the absorption. Additionally, the fact that the 
EWs of the Fe xxvi line measured here are consistent with 
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that observed with ASCA (Ueda et al. 2001) a few years 
ago, again supports this conclusion. 
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